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Scope  of  This  Talk 


>  Overview  of  Univ.  Florida  pAV 

>  Summarize  computational  capabilities  for 
fluid/structure  interactions:  membrane  and 
surrounding  viscous  flow. 

>  Present  the  aerodynamics  of  self-excited  membrane 
and  MAV  implications. 

>  Discuss  the  wing  shape  optimization  for  pAV 
Applications. 
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Gross  Weight  (Lbs.) 


Characteristics  of  jxAV 


ISi^  >  Micro  Air  Vehicle  (pAV) 

C-5 

1  smaller  than  6",  Speed  lOm/s. 
Many  applications. 

>Low  Reynolds  number  (104- 
1 05 )  condition:  degraded  L/D 

> Flight  environment 
intrinsically  unsteady. 
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The  Great  Flight  Diagram 
(modified  from  Tennekes) 
weight  ~  /3 
wing  loading  ~  / 
wing  beat  freq.  ~  /'* 
stall  speed  ~  (2W/rSCv) 

—  / 

P/W=  (D/L)V  ~  l  0  5 

Small  Birds: 

Can  fly  slower, 

Need  to  flap  faster, 

Need  less  energy  density, 
But  can  store  MUCH  less, 
Can  sustain 
higher  impact  velocity. 


jiAV :  Geometric  &  Aerodynamic 

Scaling 


Geometric  Scaling:  If  aerodynamics  is 
unchanged,  the  power  requirement 
decreases  as  the  vehicle  size  is  reduced. 

Aerodynamic  Scaling: Aerodynamic 
performance  degrades  as  the  vehicle  size, 
and  hence  Re,  decreases. 
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Low  Reynolds  Number  Airfoils 


•Gusts  affect  small  birds  and  jjAVs  more  than  large  ones 
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Representative  Low-Reynolds-Number 
Airfoils  (from  Lissaman) 


Selected  Airfoil  Profiles 
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Effects  of  Re,  Airfoil  Shape, 
and  AoA  on  Power  Index 


Membrane-Based  /iAV  Concept  at  U.  Florida 


Wingspan:  6  inches 
> Length:  5.5  inches 

>  Weight  w/payload,  video  camera: 

2  ounces 

>  Range:  0.5  mile  with  off  the  shelf 
components 

> Endurance:  10  minutes 

> Speed  Range:  10-35  miles/hour 

> Propulsion:  electric  motor 

> Batteries:  rechargeable  Lithium 
polymer 


> Altitude:  up  to  500  feet  AGL 
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Bat  Wing  Morphology 


Membrane  wing 


Adaptive  Washout  for 
Gust  Suppression 


How  direction  ^Mng  Prior 
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Computational  Fluid/Structure  Interaction 


>  Fluid  solver:  Calculate  the  external  force. 

>  Structure  solver:  Calculate  the  shape  change. 

>  Moving  boundary:  Regenerate  the  CFD  grid 

>  Interface:  Exchange  information  between  fluid  and  structure  solvers. 


Approach 

>  Structure  model 

■  Dynamic  membrane  model  with  finite  element. 

■  Expect  substantial  deformation:  nonlinearity. 

>  Fluid  flow  solver 

■  A  pressure-based  method  for  3-D  full  Navier-Stokes 
equations 

>  Grid  regeneration 

■  3 -stage  algebraic  TFI-like  method. 

>  Interpolation 

■  Thin  Plate  Spline  (TPS)  interpolation  method 
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Displacement  of  trailing  edge  at  mid-span 


Steady  Free  Stream, 
Re=  9xl04,  AoA=6° 

Periodic  oscillation  of 
the  trailing  edge  point. 

Frequency=67Hz. 
(Typical  wind  gust:  1  Hz) 

A  The  effective  angle  of 
attack  reduced. 
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Wing  Cross  Section:  Optimization? 

Pigeon  wing  Conventional  Wing 


Wing  root 


Mid-span 


Near  tip 


Optimization  Scope  and  Approach 


Minimize  CD  /  CL 

Subject  to 


1:  CT  >C 


Zbaseline 


2:  Convexity  constraint: 


Yt> 


+  .V2  “  h> 


(Xi-X0)  +  >’o-Ji+f 


x2  x0 

3:  YtL  <Yt<  i  =  l,N 

Maximize  L/D;  Maintain  lift;  Keep  cross-section  convex. 

A  direct  optimization  of  membrane  wing  is  time-demanding:  Optimize  the 
rigid  wing  as  a  surrogate. 

Design  Optimization  Tools  (DOT)  used  as  the  optimizer. 


An  automatic  grid  regeneration  tool  is  used  to  regenerate  the  CFD  grid  as 
each  analysis. 
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Choice  of  Design  Variables 


The  baseline  design  is 
based  results  from  Xfoil 
(Drela):  which  uses  a  two- 
equation  boundary  layer 
integral  formulation  & 
inviscid-BL  coupling. 

6  Design  Variables:  Three 
each  on  battens  1  and  2. 
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Airfoil  Shapes  in  Spanwise  Direction 


span  station  Z=0  span  station  Z=Q .40  span  station  Z=Q.8G 


Compared  to  the  baseline,  camber  decreases  near  the  root 
while  increases  near  the  tip. 


Overall,  the  camber  is  still  higher  at  the  root  (4.8%)  than  at 
the  tip  (4%). 


Ain  optimization  we  maintain  angle  of  attack  at  6°. 
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Spanwise  Aerodynamics  at  Design  Point: 

Rigid  Wing  at  AoA=6° 


The  improvement  is  largely  located  within  70%  of  the  inner  wing. 
r  Lift  coefficient  maintains  the  same  even  though  camber  reduces. 
The  improvement  is  largely  due  to  lower  form  drag. 


•  Department  of  Mechanical  and  Aerospace  Engineering 
•  Computational  Thermo-Fluids  Group  • 


Velocity  Profile  Near  Root:  Rigid  Wing 


>AoA=6° 

The  optimized  wing 
suppresses  the  flow 
separation. 
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Spanwise  Aerodynamics  at  Off-Design  Point: 

Rigid  Wing  at  AoA=3° 


^The  improvement  is  substantial  at  low  AoA,  and  consistent  with  the 
design  point,  is  largely  located  within  70%  of  the  inner  wing. 


Same  as  the  design  point,  the  lift  maintains  the  same  even  though 
camber  reduces,  and  the  improvement  is  largely  due  to  lower  form 
drag. 
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Aerodynamics  Between 
Membrane  &  Rigid  Wings 


Membrane  Wing  Performance  Rigid  Wing  Performance 


Optimized  shape  improves  L/D  consistently. 

Optimized  membrane  wing  varies  less  in  L/D  versus  AoA. 
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Optimized  Membrane  Wing  at  AoA=6° 


While  there  seem  substantial  variations  in  time,  the  frequency 
(about  70Hz)  is  higher  than  that  environmental  fluctuation  or  vehicle 
response. 
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Outstanding  Issues/Opportunities 


Optimized  materials  properties  for  passive 
flow  control. 

Sensor  and  simplified  aerodynamic  model 
to  facilitate  autonomous  flight  control. 

Detailed  wind  tunnel  measurements  and 
numerical  simulations  to  assess  the 
unsteady  flight  environment. 

Efficient  propulsion. 


